Accepted for publication in the Astrophysical Journal 



A NEW CEPHEID DISTANCE TO THE MASER-HOST GALAXY NGC 4258 
AND ITS IMPLICATIONS FOR THE HUBBLE CONSTANT^ 

L. M. Macri^ 

National Optical Astronomy Observatory and 
950 North Cherry Ave., Tucson, AZ 85719, USA 



K. Z. Stanek 
Department of Astronomy, The Ohio State University and 
140 West 18th Ave., Columbus, OH 43210, USA 



O 
O 

(N 
< 



(N 
00 

O 

o 
o 



D. Bersier 

Astrophysics Research Institute, Liverpool John Moores University and 
Twelve Quays House, Egerton Wharf, Birkenhead CH41 ILD, England 

L. J. Greenhill 

Harvard-Smithsonian Center for Astrophysics and 
60 Garden St., Cambridge, MA 02138, USA 

AND 

M. J. Reid 

Harvard-Smithsonian Center for Astrophysics and 
60 Garden St., Cambridge, MA 02138, USA 

ABSTRACT 

We present initial results from a time-series BVI survey of two fields in NGC 4258 using the Advanced 
Camera for Surveys onboard the Hubble Space Telescope. This galaxy was selected because of its 
accurate maser-based distance, which is anticipated to have a total uncertainty of 3%. The goal 
of the HST observations is to provide an absolute calibration of the Cepheid Distance Scale and to 
measure its dependence on chemical abundance (the so-called "metallicity effect"). 

We carried out observations of two fields at different galactocentric distances with a mean abundance 
difference of 0.5 dex. We discovered a total of 281 Cepheids with periods ranging from 4 to 45 days 
(the duration of our observing window). We determine a Cepheid distance modulus for NGC 4258 
(relative to the LMC) of A/ip =10.88 ± 0.04 (random) ±0.05 (systematic) mag. Given the published 
maser distance to the galaxy, this implies ^q{LMC) = 18.41 ± 0.10^ ± 0.13^ mag or D{LMC) = 
48.1 ± 2.3r ± 2.9s kpc. We measure a metalhcity effect of 7 = -0.29 ± 0.09^ ± 0.05^ mag dex~^ 
We see no evidence for a variation in the slope of the Period-Luminosity relation as a function of 
abundance. 

We estimate a Hubble Constant of Hq — 74 ± 3^ ± 6s km s^^ Mpc^^ using a recent sample of 4 
well-observed type la SNe and our new calibration of the Cepheid Distance Scale. It may soon be 
possible to measure the value of Hq with a total uncertainty of 5%, with consequent improvement in 
the determination of the equation of state of dark energy. 

Subject headings: Cepheids — distance scale — galaxies: individual (NGC 4258) 



1. INTRODUCTION 

During the last 15 years, the Hubble Space Telescope 
{HST) has been used to discover ~ 10'^ Cepheid variables 
in ^ 30 galaxies with D < 25 Mpc, mostly through V- 
and /-band observations carried out with the WFPC2 
instrument. The distance moduli to these galaxies have 
been determined through the use of a fiducial Cepheid 
Period-Luminosity relation (P-L) based on observations 
of variables located in the Large Magellanic Cloud. Sev- 
eral secondary distance indicators (such as type la SNe, 
the Tully-Fisher relation, the Surface Brightness Fluc- 
tuation method) have been calibrated based on these 
Cepheid distances. As a result of these investigations, 
there is some agreement that Hq is about 70 km s^^ 

^ Based on observations with the Advanced Camera for Sur- 
veys onboard the NASA/ESA Hubble Space Telescope, obtained 
at STScI, which is operated by AURA, Inc., under NASA contract 
NAS 5-26555. These observations are part of program # GO-9810. 
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Mpc""'^, perhaps with as little as 10% uncertainty (Freed- 
man et al. 2001). However, two significant sources of 
systematic error stand out. 

First, the entire Cepheid Distance Scale is underpinned 
by the distance to the Large Magellanic Cloud (LMC). 
The distance to that galaxy is used to establish the ab- 
solute calibration of the Cepheid P-L relations, and its 
uncertainty dominates the calibration of any secondary 
distance indicator. The suitability of the LMC for this 
purpose is problematic, since independent estimates of its 
distance disagree by as much as 0.5 mag, or 25% (Bene- 
dict et al. 2002). Additionally, the internal structure of 
the galaxy along the line of sight remains poorly under- 
stood (Nikolaev et al. 2004; van der Marel 2001). Faced 
with this situation, most Cepheid-based determinations 
of Hq have adopted fiLMC = 18.5 ± 0.1 mag, which cor- 
responds to a distance of Dlmc = 50.1 ± 2.3 kpc. 

Second, the effect of metal abundance on the Cepheid 
P-L relation is controversial. Several independent mcth- 
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ods for an observational determination have yielded a 
variety of results (Sasselov et al. 1997; Koclianck 1997; 
Kennicutt et al. 1998; Sakai et al. 2004) with the opposite 
sign to what has been predicted by some theoretical in- 
vestigations (Fiorentino et al. 2002), which also suggest 
a sensitivity to helium as well as metal content. Further- 
more, the use of V and / photometry alone in previous 
HST surveys makes it difficult to disentangle the effects 
of reddening and metallicity and adds uncertainty to the 
determination of Cepheid distances. 

We wish to establish a new Cepheid Distance Scale an- 
chor galaxy, NGC 4258, for which accurate geometric es- 
timates of distance are available. Herrnstein et al. (1999) 
estimated its distance modulus to be 29.29 ± 0.09^ ± 
0.12s mag, and it is anticipated that Humphreys et al., 
(in prep.) will reduce the total uncertainty of that esti- 
mate to < 3%. Our goal is even more compelling in light 
of the recent WMAP results (Spergel ct al. 2006) because 
many cosmological parameters depend sensitively on Hq 
(e.g. Eisenstein & White 2004; Tegmark et al. 2004; Hu 
2005). An accurate geometric distance to NGC 4258 can 
also be used to directly calibrate secondary distance in- 
dicators, such as the tip of the red giant branch (TRGB). 

This paper contains the first results of our project: 
deep time-series BVI photometry of two fields in 
NGC 4258 and the discovery and analysis of Cepheid 
variables. The paper is organized as follows: §2 con- 
tains details of the observations, data reduction and pho- 
tometry, and the search for variables; §3 describes the 
selection criteria and the Cepheid samples; §4 presents 
the determination of a Cepheid distance to NGC 4258, a 
measmcmcnt of the metallicity dependence of Cepheid- 
based distances, and a discussion of our results. 

Throughout the paper, we denote random (statistical) 
uncertainties with a subscript r and systematic uncer- 
tainties with a subscript s, i.e., ±0.10^ ± O.lOg mag. 

2. OBSERVATIONS, DATA REDUCTION AND 
PHOTOMETRY 

2.1. Observations 

We used the Hubble Space Telescope (GO program 
9810) to observe two fields located at widely different 
galactocentric radii within the disk of NGC 4258 (M106). 
This spiral galaxy is one of the brightest members of the 
Coma-Sculptor Cloud (TuUy & Fisher 1987), and has 
been classified as SAB(s)bc II-III (de Vaucoulcurs et al. 
1991) and Sb(s) II (Sandagc & Tannnann 1987). 

We carried out the observations using the Advanced 
Camera for Surveys / Wide Field Camera (Ford et al. 
2003), which consists of 2 back-illuminated SITe 2048 x 
4096 pixel CCDs. The average plate scale of the focal 
plane is 0'.'05 pix^^, making each image ^ 202" on the 
side. At the nominal NGC 4258 distance of ~ 7.2 Mpc, 
this translates to a physical size of ~ 1.7 pc pix^^. 

Hereafter, we refer to the two fields based on their 
galactocentric radii as "inner" and "outer" . The fields 
are centered at {a,d) = 12'^18'"47.518% +47°20'20'.'10 
(inner) and 12hl9'"23.89P, -(-47°11'37'.'61 (outer), in 
J2000.0 coordinates. Figure 1 shows the location of these 
fields in the context of a digitized POSS-H image of 
NGC42583. 

^ The Second Palomar Observatory Sky Survey (POSS-II) was 
carried out by the California Institute of Technology with funds 



The fields were observed on twelve separate epochs be- 
tween 2003 December 5 and 2004 January 19. The spac- 
ing of the visits followed a power-law distribution to min- 
imize aliasing (Madore & Freedman 2005). Table 1 con- 
tains a log of the observations. The fields were imaged in 
three colors on two consecutive orbits during each visit, 
following a standard two-point dither pattern that min- 
imizes the effects of the geometric distortion present in 
ACS while ensuring a robust cosmic-ray rejection. To- 
tal exposure times per epoch were 2 x 900s using the 
F435W filter (similar to Johnson B),2 x 800s using the 
filter (similar to Johnson V) and 2 x 400s using 
the F814VF filter (similar to Kron-Cousins /). One of 
the visits to the outer field was hampered by guide-star 
problems, reducing the total number of useful images by 
one relative to the inner field. 

2.2. Data Reduction and Photometry 

The raw observations were processed by the standard 
on-thc-fly-rcproccssing STScI ACS calibration pipeline, 
as described in the ACS Data Handbook (Pavlovsky et al. 
2005). Briefly, the pipeline performs bias level correction 
and subtraction, dark image subtraction, flat fielding cor- 
rection, and generation of ancillary data quality informa- 
tion. The calibrated images were downloaded from the 
STScI Archive and further processed using STSDAS and 
PyRAF^. Specifically, we used the PyDrizzle task to ap- 
ply the filter-dependent geometric distortion correction 
to individual images. 

We performed PSF photometry using the 
DAOPHOT/ALLSTAR/ALLFRAME suite of programs (Stetson 
1987, 1994) and following the general data reduction 
and analysis precepts of the HST Key Project on the 
Extragalactic Distance Scale (e.g.. Stetson et al. 1998). 
All the programs whose names appear in CAPITALS were 
developed and kindly provided to us by Peter Stetson. 

We defined the PSF as a quadratically-varying Moffat 
function with /3 = 1.5 and a fitting radius of 2 pixels. 
The PSF extended out to a radius of 10 pixels (0'.'5) 
and the local sky annulus was established from 20 to 25 
pixels (1 — 1'.'25). Aperture photometry was measured at 
logarithmically spaced radii from 3 to 10 pixels. 

2.2.1. Determination of template PSFs 

As expected in extragalactic Cepheid observations, our 
fields are rather crowded and lack bright, isolated stars 
suitable for the determination of the point-spread func- 
tion (PSF). Given the long-term stabiHty of HST, it is 
feasible to determine suitable PSFs from high S/N obser- 
vations of dense yet microwded stellar fields. To this end, 
we retrieved a observations of an outer field of the globu- 
lar cluster NGC 104 (program 9018), originally obtained 
for calibration purposes, from the HST Archive. 

We analyzed 18 images in each of the FASbW, F555W 
and FSliW filters, taken at various offset positions. We 
used 800 bright stars present in all the images to derive 

from NSF, NASA, the National Geographic Society, the Sloan 
Foundation, the Samuel Oschin Foundation, and the Eastman Ko- 
dak Corporation. The Oschin Schmidt Telescope is operated by 
the California Institute of Technology and Palomar Observatory. 
The Digitized Sky Survey was produced at the Space Telescope 
Science Institute under U.S. Government grant NAG W-2166. 

STSDAS and PyRAF are products of the Space Telescope Sci- 
ence Institute, which is operated by AURA for NASA 
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the PSF for each filter and to map its positional variation 
across the field of view of ACS/WFC. This was achieved 
using MULTIPSF, which is identical to the PSF routine in 
DAOPHOT but uses stars in multiple images. 

2.2.2. PSF photom,etry 

We carried out PSF photometry separately for each 
combination of field and filter, as follows: 

1. Image registration & master image: We used 
DAOPHOT to detect bright stars in the individual images 
and ALLSTAR to perform preliminary PSF photometry 
and obtain accurate coordinates of those objects. We 
used DAOMASTER to determine coordinate transformations 
for every image, taking as reference the first image of 
each set. We used MONTAGE to create a "master" image 
by interpolating the individual frames and applying a 
median filter. Figures 2 & 3 are color composites of the 
ACS fields, created from the master BVI images. 

2. Master object lists: We detected objects in the 
master image following a two-step iterative approach (de- 
tect all objects, subtract them from the image, detect all 
remaining objects and add them to the initial list). At 
each iteration, we required a 3(7 detection. The total 
number of objects detected were ^ 2,3 and 6 x 10^ in 
BVI , respectively, for the inner field and ~ 0.6, 1 and 
2 X 10^ in BVI , respectively, for the outer field. 

3. PSF Photometry: Wo used ALLFRAME to measure 
the magnitudes of every star in each master list across 
all individual images in a given field/filter combination. 
We used the same template PSFs (§2.2.1) for all frames 
in a given filter. This generated a total of ~ 2.3 x 10^ 
photometric measurements. 

4. Secondary standards: We searched the star lists for 
bright, isolated stars suitable to serve as secondary stan- 
dards. We identified ~ 30 — 50 suitable stars depending 

on the field and filter. We subtracted all other objects 
from each individual image and performed aperture pho- 
tometry on these secondary standards to generate curves 
of growth. These were analyzed using DAOGROW (Stetson 
1990) and compared to the curves of growth determined 
from bright, isolated stars in NGC 104. We found no dif- 
ference between the two sets and decided to use both in 
our analysis. 

5. Growth curves: We used COLLECT to apply the 
curves of growth and determine aperture corrections us- 
ing the secondary standards. The corrections were small, 
with average values of -h0.02±0.04, -0.04±0.03, -0.05± 
0.02 mag in BVI, respectively. Epochs #10 and #11 had 
substantially larger aperture corrections (~ 0.2 mag), 
probably due to telescope de-focusing. 

6. Zeropoints: We used CCDAVE to compute mean in- 
strumental magnitudes for the secondary standards. A 
typical r.m.s. scatter for these stars was 0.035 mag, 
and as low as 0.015 mag for the brightest objects. 
We used TRIAL to compute the final zero-point cor- 
rections for each frame, using as a reference the mean 
aperture-corrected instrumental magnitudes of the sec- 
ondary standards. 

7. Astrometry: We calculated celestial coordinates for 
all objects using the WCSTools/xy2sky program (Mink 

2002) and the astrometric solutions provided by STScI 
in the FITS headers of the first FbbbW image of each 
field. 



2.3. Photometric calibration 

As a first step in our photometric calibration, we 
corrected the instrumental magnitudes for the effect of 
charge transfer efficiency (CTE) by applying Eqn. (2) of 
the ACS Instrument Status Report 04-06 (Riess & Mack 
2004). We used the parameters derived by those authors 
for an aperture of 3 pixels in radius, since that is similar 
to the effective radius of the PSF for the bandpasses of 
interest (A. Riess, priv. comm.). Given the appreciable 
sky background in both inner and outer fields, the CTE 
correction was fairly small (~ 0.02 mag). 

We calibrated our photometry following the proce- 
dures of Appendix B of Sirianni et al. (2005), using the 
zero-points and color terms listed under the "observed" 
columns of their Table 22: 

V = F555W - ACQW + 25.704 - 0.054 {V - I) (la) 

V = F?>?>?>W - ACQbv + 25.701 - 0.056 {B - V) (lb) 
I = F814T4^- AC05/ -f 25.495 - 0.002 (F - /) (2) 
B = F435iy - AC05b + 25.842 - 0.089 {B - V) (3) 

where BVI are the standard magnitudes and FAZfiW, 
FbbbW, ^"81414^ arc the CTE-corrected, aperture- 
corrected (to 0'.'5) instrumental magnitudes derived in 
§2.2.2. The additional filter-dependent aperture correc- 
tions to infinity, ACQbi, are listed in Table 5 of Sirianni 
et al. (2005). In the case of objects with three-color pho- 
tometry (such as all the Cepheid variables) , we gave pref- 
erence to (la) over (lb) because our /-band observations 
have higher S/N than our B-band observations. Since 
these transformation equations make use of the standard 
(rather than observed) colors, they were applied itera- 
tively until convergence. 

Table 2 lists the positions and calibrated magnitudes of 
the secondary standards to facilitate future comparisons 
with our work. 

2.4. Search for Variables and Classification 

We searched for variables using the TRIAL program, 
which performs a scaling of the reported ALLFRAME mea- 
surement errors and calculates robust mean magnitudes 
and modified Welch-Stetson variability indices Ly (Stet- 
son 1996). Figure 4 shows the distribution of Ly as a 
function of V magnitude for the outer field (the inner 
field distribution is very similar but denser). By con- 
struction, the mean value of Ly is zero. Given the ob- 
served Icr dispersion in Ly of 0.25, we set L > 0.75 as 
the minimum variability threshold. We calculated the 
twenty most likely periods for each variable using the 
Lafler-Kinman algorithm (Lafler & Kinman 1965) as en- 
coded in TRIAL. 

We applied an automated classification algorithm de- 
veloped by the DIRECT project (Kaluzny et al. 1998) 
to the F-band light curves of the variables. The algo- 
rithm computed the chi-squared per degree of freedom, 
X^, of each light curve for three cases: i) a constant mag- 
nitude (null hypothesis) ; ii) a linearly- varying magnitude 
(appropriate for objects with periods much longc;r than 
our observing window); iii) a Cepheid variable with a 
period equal to each one of the twenty tentative periods 
returned by the Lafler-Kinman algorithm. The latter 
case used the Cepheid template light curves developed 
by Stetson (1996). 
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We selected as possible Cepheids those periodic vari- 
ables with a xi foi' case (iii) that was at least a factor 
of two smaller than the xi of cases (i) or (ii), follow- 
ing the methodology of the DIRECT project. Adopting 
the best-fit period, we phased the B- and /-band light 
curves and fit them with Cepheid template light curves 
(absent _B-band templates, we used suitably scaled V- 
band ones). Finally we phase- weighted mean magnitudes 
(Stetson 1996) through numerical integration of the best- 
fit template light curve for each variable. 

3. SELECTION OF CEPHEID VARIABLES 

3.1. Adopted Period-Luminosity relations 

Our analysis requires the adoption of fiducial Period- 
Luminosity (P-L) relations to calculate distance moduli, 
as well as corresponding Period-Color (P-C) relations 
to correct for the effects of interstellar extinction. We 
adopted the P-L relations originally derived by Udalski 
et al. (1999), as updated in the OGLE web site'"'. These 
relations are based on a sample of iV > 600 Cepheids 
observed as part of the OGLE II project, with periods 
ranging from 2 to 30 days: 

B = 14.929(31) - 2.439(46) [log P - 1] (4) 
V = 14.287(21) - 2.779(31) [log P-1] (5) 
I = 13.615(14) - 2.979(21) [log P-1] (6) 

where P is the period of the Cepheid in days, and the 
errors in the zero-points and slopes are expressed in units 
of 10"'^. The dispersions of the data relative to the rela- 
tions are 0.24, 0.16 and 0.11 mag in BVI, respectively. 

3.2. Extinction corrections and relative distance moduli 

NGC4258 is located at I = 138°32, b = 68°84. We 
estimated the value of foreground Galactic interstellar 

extinction based on the values in the reddening map of 
Schlegel et al. (1998) for a number of positions near the 
galaxy. All of them yielded very low values of foreground 
extinction, E{B — V) = 0.016 mag. We expect little 
additional (internal) extinction in the outer field, but 
the Cepheids in the inner field should be subject to a 
considerably larger amount of internal extinction with 
strong variations as a function of position. 

We determined the total extinction to each Cepheid 
by comparing the observed B — V, V—I and B — I col- 
ors with the zero-extinction colors {B — V)q, {V — I)o 
and (B — I)q predicted by the P-C relations formed by 
Eqns. (4-6). We transformed the values of E{V — I) and 
E{B-I) to E{B-V) using the values of from Table 6 
of Schlegel et al. (1998) for Ry = 3.1 and the extinction 
law of Cardelli et al. (1989). We used the three color 
excesses to compute a mean E{B — V) and standard de- 
viation, but adopted a 0.025 mag uncertainty floor to 
accoimt for the intrinsic width of the P-C relations. 

We determined extinction-corrected LMC-relative dis- 
tance moduli for each Cepheid by calculating: 

A/^o = Am/ - IA^E{V-I) (7) 
Afj.a^Ani-2ME{B-I) (8) 
Afio==Afxi-lME{B-V) (9) 

^ ftp;//sirius. astrouw.edu.pl/ogle/ogle2/vaj_stars/lmc/ 
cep/catalog/README.PL 



where A/j,/ is obtained by subtracting the mean value 
of I{P) from Eq. (6) from the mean /-band magnitude 
of a given Cepheid of period P. The values of total- 
to-selective extinction ratios were calculated using the 
Ax/E{B-V) values in Table 6 of Schlegel et al. (1998). 
The three values of A/iq were averaged to obtain a mean 
value and standard deviation. Eqn. (7) is the classical 
Wesenheit distance modulus used by Freedman et al. 
(2001) and Saha et al. (2001), although these authors 
apply it to distance moduli determined from an ensem- 
ble of Cepheids. We chose to form the additional two 
relations (Eqns. 8-9) to rule out any global systematic 
error in the calibration of the photometry. 

We note that in this analysis, all the calculated dis- 
tance moduli are relative to the LMC, since the adopted 
P-L relations are expressed in observed magnitudes. 
We adopt this approach to accommodate the antici- 
pated improvement in the geometric distance estimate 
for NGC4258 to be obtained from new analyses of the 
maser observations (Humphreys et al., in prep.). 

3.3. Selection criteria 

The selection of Cepheids from a larger set of periodic 
variables is not a trivial undertaking for HST observa- 
tions, especially given the crowded nature of the fields, 
the sparse sampling of the light curves and the relatively 
low S/N of the individual data points at the faintest 
magnitudes (corresponding to the shortest periods). Dif- 
ferent studies have adopted various selection techniques, 
some based on visual inspections of light curves and im- 
ages (Saha et al. 1996) and some based on a more math- 
ematical approach (Leonard et al. 2003). 

In the case of these observations of NGC 4258 — which 
is located significantly closer than most Cepheid-bearing 
galaxies studied with HST~ it is fairly easy to select a 
sample of high-quality Cepheids with P > 10 days for 
the primary scientific goals of deriving a distance and 
measuring the metallicity effect. The values we obtain 
are insensitive to the application of different selection 
techniques. The selection of shorter-period (P < 10 d) 
Cepheids is less certain, especially in the case of the in- 
ner field. We have adopted a particular set of selection 
criteria, but we list all detected Cepheid candidates to 
facilitate alternative analyses by others. 

We restricted the sample of variables using two cuts in 
Lv' a fairly low value of 0.75 (hereafter, the "extended 
sample" ) and a more conservative value of Ly = 2 (here- 
after, the "restricted sample"). We applied additional se- 
lection criteria based on observed properties (light curve 
amplitude ratios, colors, etc.) in an attempt to remove 
contaminated Cepheids from the sample. 

We detected a total of ~ 10^ distinct objects at the 
> 3(7 level in the master images. To exclude false pos- 
itives, we trimmed objects detected in less than 75% of 
the individual images of each band. 

As stated in §2.4, variable stars were classified as 
Cepheids if the xt '^a-s reduced by more than a factor 
of two (relative to the null hypothesis) when fitting the 
phased T^-band data with a template Cepheid light curve. 
These requirements were met by 536 objects in the ex- 
tended sample and 258 objects in the restricted sample. 
We applied the following secondary selection criteria: 
1. Amplitude ratios: Relative BVI amplitudes for 
fundamental- mode pulsators obey the proportions 1.5 : 
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1 : 0.5. To discard obvious blue or red blends, we re- 
quired (a) an I-to-V amplitude ratio of 0.5 ± 0.25 and 
(b) a B-to-V amplitude ratio of 1.5 ± 0.5. 

2. Blue edge: We discarded objects with E{B — V) 
more than 2a below the Galactic foreground value of 
0.016 mag, as these variables are likely to be blended 
with blue stars. Wc set E{B — V) = 0.016 mag for ob- 
jects between the threshold and the foreground value. 

3. Large extinction: We rejected objects with 
E{B — V) > 0.5 mag. These Cepheids are either blended 
with red stars or highly reddened, in which case the ac- 
tual value of Rv could deviate significantly from the 
adopted value of 3.1. 

4. Pop II Cepheids: We imposed a conservative up- 
per limit of A/io < 12 mag to the LMC-relative distance 
moduli to reject long-period RV Tauri and W Virginis 

variables in our sample. These Population II Cepheid- 
like variables obey period-color relations that are similar 
to those of Cepheids, and therefore had passed the pre- 
vious selection criteria. However, their absolute magni- 
tudes are more than 1 mag fainter than Cepheids. 

Once these cuts were applied, we computed a mean 
modulus for each sample using a least-absolute-deviation 
technique with iterative sigma clipping. This was moti- 
vated by the asymmetric tails of outliers that are caused 
by other sources of contamination, many of which make 
objects artificially brighter. 

3.4. Final Cepheid samples 

Table 3 details the effects of the selection criteria (§3.3) 
on the initial samples. 281 out of 536 variables in the 
extended sample pass all criteria, while 173 out of 258 
variables in the restricted sample remain. Figure 5 shows 
the period distributions for both cuts in Ly- We used the 
restricted sample (Ly > 2) in the subsequent analysis. 

The observed properties of the 281 Cepheids that 
passed our selection criteria are listed in Table 4, while 
those of the 255 rejected candidates are listed in Table 5 
along with the reasons for their rejection. Table 6 lists 
the derived properties of the Cepheids in Table 4. Ta- 
ble 9 contains the individual photometric measurements 
of these objects. 

Figures 6 & 7 show the distribution of the Cepheids 
within the outer and inner fields, respectively; individual 
finding charts can be see in Figures 8a-g. Figures 9 & 10 
show the distribution of the Cepheids within the color- 
magnitude diagrams of the two fields. Representative 
light curves are shown in Figures 11 & 12. Lastly, Fig- 
ures 13 and 14 contain the observed BVI P-L relations 
of the restricted samples. 

4. RESULTS AND DISCUSSION 
4.1. The Maser Distance to NGC 4258 

Water maser emission observed from NGC 4258 orig- 
inates in a subparsec annular region within a nearly 
edge-on, warped accretion disk, bound by a supermassive 
black hole in the nucleus (Miyoshi et al. 1995; Greenhill 
et al. 1995). Masers lie: (1) in a narrow sector on the 
near side of the disk and (2) on the disk-diameter per- 
pendicular to the line of sight. 

Geometric estimates of distance may be obtained 
from measurements of the centripetal acceleration or the 
proper motion of masers on the near side of the disk. 



The acceleration is obtained from the time rate of change 
of the maser Doppler shifts, and the proper motion is 
obtained from the change in the positions of the near- 
side masers relative to the approaching/ receding masers 
(which appear to be stationary on the sky). 

Herrnstein et al. (1999) reported acceleration and 
proper motion distance moduli that agreed to < 1%: 
Mmaser = 29.29 ± 0.09^ ± 0.12^ mag. The quoted system- 
atic uncertainty arises largely from unmodeled structure 
and an upper limit on the eccentricity of the disk. Initial 
models assumed circular orbits and a warp in position 
angle alone. More recently, Herrnstein et al. (2005) per- 
formed a detailed analysis of the disk rotation curve and 
detected a 2a deviation from a Keplerian law, which they 
attributed to an inclination-warp in the disk. That also 
helps to explain the locus of the near-side masers. 

Humphreys et al., (in prep.) aim to reduce the random 
component of the uncertainty by including more epochs 
of observation, and more importantly, to reduce the sys- 
tematic component by improving the dynamical model 
of the maser-disk system. The Herrnstein et al. (1999) 
distance relied on VLBI data collected at four epochs 
between 1994 and 1997, while data for 18 VLBI epochs 
(1997-2000) and 40 spectroscopic epochs (1994-2003) are 
now available. The analysis also limited disk eccentric- 
ity to < 0.1. More densely sampled data with a longer 
time baseline, coiipled with a more sophisticated model 
of the disk warp and eccentricity, are anticipated to re- 
duce the systematic and random uncertainties in distance 
by more than a factor of two, for a total uncertainty of 
- 3% (Humphreys et al. 2005a,b). 

4.2. A Cepheid distance to NCC4258 

4.2.1. Minimum period cut 

We imposed minimum period cuts to the samples de- 
rived in §3.4 before we determined mean relative distance 
moduli. Several reasons motivate the use of such a cut. 

a) We are unable to differentiate between fundamental 
and overtone pulsators due to our sparse phase sam- 
pling. Overtone pulsators in the Magellanic Clouds have 
2 < P < 6 days and are ~ 0.75 mag brighter than fun- 
damental pulsators with the same period (Udalski et al. 
1999). Hence, they can produce a large systematic bias 
in the derived distance. 

b) Confusion noise introduces a systematic bias in the 
photometry of Cepheids that becomes increasingly im- 
portant at faint magnitudes, especially in the /-band 
(Saha & Hoessel 1990; Saha et al. 1996). 

c) Observing objects near the detection limit may result 
in incompleteness bias at the shortest periods of the ob- 
served P-L relation (Sandage 1988). 

d) The observed magnitudes of short-period Cepheids 
are more likely to be contaminated by unresolved blends 
with other disk stars (Mochejska et al. 2000), especially 
in the denser regions of the inner field. 

We applied the cut at minimum period and calculated 
the mean value of the individual relative distance moduli 
following the procedure described in §3.3. Figures 15 
& 16 show the impact of this procedure for the outer and 
inner fields. Figure 17 shows the mean relative distance 
modulus and its uncertainty as a function of Pmin for 
both fields. 

There is no statistically-significant variation in the 
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mean relative distance modulus of the outer field as a 
function of Pmin- The primary use of the outer field 
Cepheids in this study is to test the maser distance of 
NGC 4258 against the distance to the LMC without hav- 
ing to worry about abundance differences (since they 
have the same mean metallicity). We chose Pmin = 6 d 
as the final period cut for this sample to avoid contami- 
nation by overtone pulsators and to maximize the sample 
size and the overlap of period ranges between these two 
galaxies; note that Pmax is 32 d for the OGLE LMC 
sample and 44 d for NGC 4258. 

The inner field exhibits a mild trend with shorter dis- 
tance moduli for smaller minimum period cut-offs, with 
a statistical significance of ^ 2.7 cr [Pmin of 20 d vs. 
5 d). We chose Pmin = 12 d to avoid the observed bias 
in distance modulus at shorter periods. This value of 
Pyyiin IS similar to the typical lower limit of the Ccpheid 
samples discovered in other galaxies observed with HST 
(P„i„ = 10-15 d). 

4.2.2. Distance moduli 

Taking the aforementioned period cuts into account, 
and using the restricted samples, we derive distance mod- 
uH relative to the LMC of A/xq = 10.87±0.05,.±0.05s mag 
(outer field, iV = 20 Cepheids) and A^o = 10.71±0.04^± 
0.05s niag (inner field. A'' = 69 Cepheids). 

The quoted uncertainties for these relative distance 
moduli arise from terms B & C of our error budget, which 
is listed in detail in Table 7. For comparison, we also list 
the error budget typical of Cepheid distance determina- 
tions based on HST/WFPC2 observations (e.g., Gibson 
et al. 2000) as well as the anticipated error budget after 
our follow-up NICMOS and ACS/HRC data are incor- 
porated in the analysis and the uncertainty in the maser 
distance is reduced. 

We derived relative distance moduli for the two fields 
using the methodology of Freedman et al. (2001), in 
which one calculates mean V and / distance moduli for 
the Cepheid ensemble (i.e., neglecting differential red- 
dening among Cepheids). We did not apply rejection 
criteria 2-|-3, since they were not used by those authors, 
and used the same period cuts as above. We obtained 
A/xy = 11.19±0.04^±0.02s, A^/ = 11.05±0.03r±0.02s 
and A/io = 10.86 ± 0.04,. ± 0.05^ mag [N = 38, outer), 
andA/iy = 11.33±0.04r, A/z/ = 11.07±0.03r and A^o = 
10.69 ± 0.04,. ± 0.05s mag [N = 85, inner). These val- 
ues are consistent with a previous HST/WFPC2 Cepheid 
distance to NGC 4258 derived by Newman et al. (2001) 
using the same methodology. Their ALLFRAME pho- 
tometry of A'' = 7 Cepheids with P = 10 - 21 d yielded 
A/xo = 10.90 ± O.lOr ± 0.06s mag. 

4.3. Metallicity dependence 

The two fields under study provide an excellent op- 
portunity to obtain a differential measurement of the 
metallicity dependence of the Cepheid P-L relation. We 
adopted an abundance gradient for NGC 4258 measured 
by Zaritsky et al. (1994) and expressed in their "empiri- 
cal" metallicity scale as: 

[O/H] = 8.97 ± 0.06 - 0.49 ± 0.08(p - 0.4) dcx (10) 

where p is the deprojected galactocentric radius, ex- 
pressed as a fraction of the isophotal radius pq: 

x = {a — ao)cos^+ {6 — So)sm^ (11) 



y = {{6 - So) cos (f)- {a- ao) sin 4>}/{b/a) (12) 

p={x'+y^)i/po (13) 

We computed the deprojected galactocentric distances 
of the Cepheids using these equations. We adopted <p = 
U9°75, hi a = 0.413 and po = 7'.76 (derived from a least- 
squares fit to the data in Table 2 of Zaritsky et al. 1994) 
and a position for the center of NGC 4258 in J2000.0 
coordinates of {a,d) = 12'^18'^57.«5046, -|-47°18'14'.'303 
(Herrnstein et al. 2005). 

Figure 18 shows the correlation between true distance 
modulus and deprojected galactocentric distance, or its 
corresponding abundance according to Eqn. 10. The 
sample plotted in this figure comprises all Cepheids from 
Table 4 with Ly > 2 (i.e., the restricted sample) and 
P > 6 d (outer field) or P > 12 d (inner field). At 
the suggestion of the referee, we further restricted the 
samples to ensure that they cover the same range of ex- 
tinction, 0.05 < E{B-V) < 0.28 mag (N=69 Cepheids). 

A least-squares fit to the data yields 7 — —0.29 ± 
0.09r. ± 0.05s mag dex and A/xo(NGC4258- 
LMC)= 10.88 ± 0.04,, ± 0.05s mag, measured at 
12 -|- \og[0/H] = 8.5 dex. The best fit is represented 
by a solid line in Fig. 18. Since this is a differential mea- 
surement within a single galaxy, the random uncertainty 
arises from the scatter in the individual distance mod- 
uli and the systematic error is due to the uncertainty in 
the determination of the Zaritsky et al. gradient. Fig- 
ure 19 shows the residuals of the individual distance mod- 
uli about the fit, plotted as a function of E{B — V). 

Our measurement compares favorably with the recent 
determination of Sakai et al. (2004), who used the tip 
of the red giant branch as a fiducial distance indicator 
under the assumption that is unaffected by abundance 
differences. They derived 7 = —0.25 ± 0.09 mag dex~^ 
by comparing distances determined using Cepheid vari- 
ables and the Tip of the Red Giant Branch (hereafter 
TRGB) to 17 nearby galaxies. The Cepheid distances 
were calculated using the same P-L relations we adopted 
(Eqns. 5 & 6). 

Our result is also consistent with, but more statisti- 
cally significant than an earlier differential determina- 
tion of the metallicity effect by Kennicutt et al. (1998), 
who found 7 = —0.24 ± 0.16 mag dex^^ based on HST 
observations of Cepheids in two fields within MlOl. 

Likewise, our findings are in agreement with the values 
of metallicity dependence derived by Kochanek (1997) 
through an analysis of Cepheid magnitudes and col- 
ors in multiple galaxies, and by Sasselov et al. (1997) 
from a differential comparison of Large and Small Mag- 
ellanic Cloud Cepheids. We find a difference in dis- 
tance modulus between the inner and outer fields of 
6fj,o = —0.15 ± 0.04 mag for a mean abundance dif- 
ference of AZ = 0.45 dex; the aforementioned stud- 
ies would have predicted A/ig = —0.15 ± 0.06 mag and 
—0.18 ± 0.08 mag, respectively. 

Adopting the Tg metallicity scale of Kennicutt et al. 
(2003), the coefficient of the metallicity dependence be- 
comes 7 = —0.49 ± 0.15r mag dex~^. 

4.4. A Tip of the Red Giant Branch distance to 

NGC'4258 

At the suggestion of the referee, we determined a dis- 
tance to NGC 4258 using the TRGB method (Lee et al. 
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1993; Sakai et al. 2004). The /-band master image of the 
outer field reaches a depth of / ^ 27 mag, which is signif- 
icantly deeper than the expected TRGB magnitude. The 
F-band master image reaches a depth of F ^ 28 mag, 
which is sufficient to reject all stars in the I-band lumi- 
nosity function with V — I<1 mag. Such a color cut is 
standard practice in TRGB studies Sakai et al. (2004); 
Mendez et al. (2002). 

The outer field /-band master object Hst (§2.2.2, 2), 
contains 2.05 x 10^ objects. We rejected objects that ap- 
peared in less than half of the individual frames or exhib- 
ited signs of variabihty (L/ > 0.75), reducing the sample 
to 1.37 X 10"'' objects. Then, we rejected a small fraction 
(2%) of the remaining objects which exhibited a poor 
fit to a stellar PSF relative to other objects of the same 
magnitude. These are either faint galaxies or crowded 
stars. Wc matched the remaining 1.35 x 10"'' objects 
against the V-band master list and rejected all objects 
with V — I<1 mag. Thus, the final /-band luminosity 
function that served as input for the TRGB detection al- 
gorithm consisted of 1.2 x 10^ stars with V—I>1 mag. 

We computed the TRGB magnitude following the pro- 
cedures described in Sakai et al. (1996) and Mendez et al. 
(2002). We computed a continuous luminosity function 
<j){m) using Equation (Al) of Sakai et al. (1996) and a 
logarithmic edge-detection function E{m) using Equa- 
tion (4) of Mendez et al. (2002). We measured the TRGB 
magnitude by identifying the highest peak in the prod- 
uct E{m)\/ (l>{m) and fitting a cubic spline to the region 
±0.15 mag about the peak. Lastly, we estimated the un- 
certainty in our measurement of the TRGB magnitude 
by performing a bootstrap test with 500 simulations, as 
carried out by Sakai et al. (2004). 

The right panel of Figure 20 shows the values of 4'{m) 
and E{m) that we obtained, resulting in a clear detection 
of the TRGB at Itrgb = 25.42 ± 0.02 mag. For refer- 
ence, the TRGB magnitude is also shown as a dashed line 
in the /-band CMD plotted in the left panel of Figure 20; 
note that the actual dataset used to measure the TRGB 
was far more complete than what can be shown in the 
CMD, containing 4x more stars with V — I>1 mag and 
reaching / ^ 27 mag. 

We corrected the observed / magnitude of the TRGB 
for foreground reddening (§3.2) by Aj — 0.03 mag. We 
also applied bolomctric and mctallicity corrections, fol- 
lowing Equations (l)-(4) of Sakai, Zaritsky & Kennicutt 
(1999). These equations require the determination of 
the mean V — I color of stars at the TRGB edge and 
0.5 mag below it. We determined those values to be 
{V-I)trgb = 2 ±0.25 and (y-/)_3.5 = 1.75 ±0.25 by 
constructing histograms of the V—I color distribution 
for stars within ±0.1 mag of / = 25.42 and 25.92 mag, 
respectively. The bolometric and mctallicity correction 
amounts to +0.02 ± 0.08 mag 

After these corrections, we find Itrgb — 25.41 ± 
0.04r± 0.08s mag. The corresponding value for the LMC 
(Sakai, Zaritsky & Kennicutt 1999) is /^^^^(/.MC) 
14.54 ± 0.04^ ± 0.06s mag. Thus, we determine an LMC- 
relative distance modulus to NGC 4258 based on the 
TRGB method, of A/xo,thgb = 10.87±0.06^±0.10s mag, 
in excellent agreement with the Cepheid relative distance 
modulus obtained in §4.3. Additionally, this determina- 
tion allows us to increase the sample of galaxy fields used 



by Sakai et al. (2004) to determine the Cepheid mctal- 
licity dependence based on the observed difference be- 
tween TRGB and Cepheid distance moduli. Figure 21 
is an updated version of the bottom panel of Figure 
12 of Sakai et al. (2004), with the addition of the two 
fields in NGC 4258. The best-fit line to the data is 
7 = —0.27 ± 0.06 mag dex^^, in very good agreement 
with the metallicity dependence we independently de- 
rived in §4.3. 

4.5. Other Period-Luminosity relations 

We considered in our analysis a second set of LMC 
Period-Luminosity relations derived by Sandage et al. 
(2004) using the sample of Udalski et al. and additional 
long-period Cepheids (P = 10 — 80 d) from the literature. 
These P-L relations have two slopes, with the break point 
set at P = 10 d (as motivated by Kanbur & Ngeow 2004). 
We found no statistically significant difference between 
the distance moduli derived using the Udalski et al. and 
the Sandage et al. relations. This is consistent with 
the observation of Ngeow et al. (2005) that very large 
samples of Cepheids {N > 10^) are required to detect 
the change in slope of the P-L relations. 

Additionally, we considered the P-L relations derived 
by Tammann et al. (2003) for Milky Way Cepheids. In 
that study, the individual distance to each variable was 
derived using the Baadc-Wcssclink method and/or the 
open-cluster main-sequence fitting method. The authors 
determined P-L relations with slopes that were signifi- 
cantly steeper than those derived using LMC Cepheids. 
They attributed the change in slope to abundance dif- 
ferences, since the Milky Way Cepheids in their sample 
have a mean metallicity that is close to solar. Recently, 
Saha et al. (2006) recalibrated the peak luminosities of 
type la SNe using P-L relations whose slopes vary as a 
function of abundance and Sandage et al. (2006) used 
the results to derive Hq = 62 ± 6 km Mpc^^. 

However, there is some controversy over the P-L re- 
lation slopes that arc derived via the Baade-Wesselink 
method. Gieren et al. (2005) applied this technique to 
LMC Cepheids and derived diff'erent P-L slopes than 
those of Udalski et al. (1999). They attributed the differ- 
ence to a systematic error in the Baade-Wesselink tech- 
nique, which requires the use of a period-dependent pro- 
jection factor p. Gieren et al. proposed a new p factor 
that would resolve the discrepancy. However it still re- 
mains to be explained why Tammann et al. derived es- 
sentially identical Milky Way P-L relations using a com- 
pletely independent method (open cluster main-sequence 
fitting). Parallax measurements to Galactic Cepheids be 
provided by GAIA in the next decade may yield a defini- 
tive answer on this matter. 

We can test the hypothesis of Saha et al. (2006) with 
our large sample of Cepheids in the inner field, since the 
application of the correct V and / P-L relations should 
lead to a distribution of distance moduli that is uncorre- 
lated with period. We started with the restricted sample 
of 195 Cepheids in the inner field and excluded 23 objects 
with anomalous amplitude ratios (§3.3.1), 24 variables 
with P < 6 d, and 5 objects with outside 11 ± 1 mag. 
Next, we fit a slope to A/xq vs P using an iterative least- 
absolute-deviation procedure with 3a clipping, which re- 
jected 6 outhers. Thus, our final sample consisted of 137 
Cepheids. We carried out this exercise for three choices of 
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P-L relation: Udalski et al. (1999), Sandage et al. (2004), 
and Tammann et al. (2003). We tested the null hypoth- 
esis by computing the Spearman rank-order correlation 
coefficient Vs for each choice of P-L relation. For compar- 
ison, we carried out the same exercise for the outer field 
sample and the P-L relations of Udalski et al. Figure 22 
shows the result of these tests. 

The LMC P-L relations are a good fit to the samples of 
both fields. There is a small correlation for the inner field 
with = 0.2 (2.5cr), which decreases to = 0.15 (1.3ct) 
if we use Pmin = 12 d as in §4.2. The application of the 
Milky- Way P-L relations of Tammann ct al. (2003) to 
the inner field sample yields a distribution that deviates 
noticeably from the null hypothesis, with = 0.6 (6.8cr). 
The correlation is still present, with Vg = 0.45 (4cr), for 
Pmin = 12 d. Thus, we conclude that the LMC P-L 
relations arc a bettor fit to both samples, regardless of 
their abundance difference. 

4.6. Implications for Hq and w 

Since the mean abundance of LMC Cepheids (12 -|- 

\og[0/H] = 8.5 dcx) lies within the range spanned 
by our sample of variables (Fig. 18) we have mea- 
sured A/xo(NGC4258-LMC)= 10.88 ± 0.04^ ± 
0.05s mag (§4.3). Combined with the maser distance 
modulus to NGC 4258, we infer the distance modulus of 
the LMC to be fio{LMC) = 18.41 ± 0.10^ ± 0.13^ mag. 
This corresponds to a distance of D{LMC) — 48.1 ± 
2.3r- ± 2.9s kpc, which is in excellent agreement with the 
value of 48.3 ±1.4 kpc derived from eclipsing binaries 
(see Case II in Table 8 of Fitzpatrick et al. 2003). Im- 
portantly, both distance estimates are mainly geometric, 
independent of each other, and do not rely on any "stan- 
dard candles" . 

We note that in the near future, there will be four 
galaxies with "geometric distances" that can serve as 
absolute calibrators for the Cepheid Distance Scale: the 
Large Magellanic Cloud (with multiple DEB distances, 
see Fitzpatrick et al. 2003, and references therein). 
Messier 31 (with a DEB distance by Ribas et al. 2005), 
Messier 33 (with a DEB distance by Bonanos et al. 2006) 
and NGC 4258 (with the maser distance by Humphreys 
et al., in prop.). Thus, we can expect a significant re- 
duction in the uncertainty of the "first rung" of the Ex- 
tragalactic Distance Scale, which has been a dominant 
source of uncertainty in recent determinations of Hq- 

The implied decrease in the distance to the LMC de- 
rived in this paper, relative to the adopted value of 
D = 50.1 ± 2.3 kpc (Freedman et al. 2001; Saha et al. 
2001), affects previously-derived values of Hq by ~ +3%. 
The increase in the coefficient of the metallicity depen- 
dence from 7 = — 0.2±0.2 mag dex~^ (adopted by Freed- 
man et al. 2001) to -0.29±0.09r±0.055 mag dex'^ (§4.3) 
has an opposite effect on Hq of ^ —2%. As a result, the 
net effect on the calibration of secondary distance indi- 
cators is mitigated. Table 8 shows a rc-calculation of the 
peak absolute V magnitude of type la SNe recently de- 
termined by Riess et al. (2005), which changes only by 
-0.03 mag to = -19.14 ± 0.06 mag. The resulting 
value of Ho is 74 ± 3^ ± 6s km 8"^ Mpc~^ 

Recently, Spergel et al. (2006) presented a determi- 
nation of cosmological parameters based on 3 years of 
WMAP observations. CMB observations cannot provide 
strong constraints on the value of Hq on their own, due to 



degeneracies in parameter space (Tegmark et al. 2004). 
Figure 23 shows the degeneracy in the Qm — w; plane. 
The addition of an independent of Hq from Cepheids 
significantly reduces that degeneracy (Hu 2005). 

We calculated the improvement due to a prior on Hq 
(solid contours of Fig. 23) by resampling the Monte Carlo 
Markov Chains kindly made available by the WMAP 
team, using Eq. B4 of Lewis & Bridle (2002). We also 
calculated marginalized probability distributions for w 
for increasingly more accurate priors on Hq. The results, 
which are shown in Figure 24, indicate that a 5% prior 
on Hq would reduce the la uncertainty in w to ±0.1. As 
shown by Spergel et al., the combination of CMB data 
with more than one prior (e.g., Cepheids, type la SNe 
and large-scale structure) can further refine the measure- 
ment of IV. 

A determination of Hq to 5% (sec Table 7) is a con- 
servative goal for the near term. It will require the re- 
estimation of a maser distance to NGC 4258 (Humphreys 
et al., in prep.), the analysis of follow-up observations of 
the Cepheids discovered in this paper with other HST in- 
struments (Bersier et al., in prep.; Macri et al., in prep.), 
and the inclusion in the Cepheid sample of longer-period 
(40 d < P <90 d) variables discovered with GMOS on 
Gemini North (Macri & Smith, in prep.). 

Further improvement on the accuracy of Hq, down to 
1%, may be obtained through maser distances to a large 
number of galaxies in the Hubble flow, which could be 
discovered with the Square Kilometer Array and its pro- 
totypes (Greenhill 2004). 

5. CONCLUSIONS 

The five main results presented in this paper are the 
following: 

1. We discovered 281 Cepheid variables in two fields lo- 
cated within the galaxy NGC 4258, with accurately cali- 
brated BVI photometry in twelve epochs per band. 

2. We determined a relative distance modulus between 
NGC 4258 and the Largo Magellanic Cloud, based on 
Cepheid variables, of A^o = 10.88 ± 0.04^ ± 0.05^ mag. 

3. We determined a relative distance modulus between 
these two galaxies, based on the Tip of the Red Giant 
Branch method, of A/^o, tkgb = 10.87 ± 0.04^ mag. 

4. We measured a metallicity dependence of the Cepheid 
distance scale of 7 = —0.29 ± 0.09^ ± 0.05^ mag dex^^. 

5. Our observations are best fit with P-L relations that do 
not invoke changes in slope as a function of abundance. 
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Fig. 1 .— a blue DPOSS image of NGC 4258 showing the two ACS/WFC fields observed for this project and the WFPC2 field previously 
studied by Newman et al. (2001). 




Fig. 2.— A color composite of the HST ACS/WFC outer field of NGC 4258. 




Fig. 3.— A color composite of the HST ACS/WFC inner field of NGC 4258. 
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Fig. 4. — Distribution of the modified Welch-Stetson variability index Ly (Stetson 1996) as a function of V magnitude for objects in the 
outer field, after sigma-scaliiig. The dotted and dashed lines represent the minimum values of Ly for the extended and restricted samples, 
respectively. Two variables with Ly > 6 are represented by arrows. 
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Fig. 5. — Period histograms for the final Cepheid samples after application of the selection criteria listed in §3.3. Left panel: Inner field. 
Right panel: Outer field. Dashed line: extended sample. Solid line: restricted sample. 




Fig. 6. — Master F-band image of the ACS outer field. The location of the Cepheids listed in Table 4 are marked with open circles. 




Fig. 7. — Master F-band image of the ACS inner field. The location of the Cepheids listed in Table 4 are marked with open circles. 
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Fig. 8b. — Individual finding charts for the Cepheids fisted in Table 4. Each box is 2.5" on a side. 
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Fig. 8g. — Individual finding charts for the Cepheids listed in Table 4. Each box is 2.5" on a side. 
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Fig. 9. — Color-magnitude diagrams for objects in the outer field of NGC 4258. Cepheids are plotted using filled and open symbols for 
the restricted and extended samples, respectively. Field stars are represented by small dots. The dashed lines represent the zero-reddening 
instability strip of LMC Cepheids and its 2(7 width. The arrows indicate the effect of E{B — V) = 0.2 mag. 
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Fig. 10. — Color-magnitude diagrams for objects in the inner field of NGC4258. Ccpheids arc plotted using filled and open circles for 
the restricted and extended samples, respectively. Field stars are represented by small dots. The dashed lines represent the zero-reddening 
instability strip of LMC Cepheids and its 2cr width. The arrows indicate the effect of E{B — V) = 0.2 mag. 




Fig. 11. — Representative light curves of Cepheids in the outer field. Blue: B\ green: V\ red: /. The sohd fines indicate the best-fit fight 
curve template from Stetson (1996). 
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Phase 

Fig. 12. — Representative light curves of Cepheids in the inner field. Blue: B; green: V; red: /. The solid lines indicate the best-fit light 
curve template from Stetson (1996). 
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Fig. 13. — BVI and Wcsenhcit Period-Luminosity relations for the Ly > 2 sample of Ccpheids in the outer field. The solid lines 
represent the LMC P-L relations derived by Udalski et al. (1999), shifted to the appropriate mean relative distance modulus. The dashed 
lines indicate the 1<t dispersion of the sample. 
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Fig. 14. — BVI and Wesenheit Period-Luminosity relations for the Ly > 2 sample of Cepheids in the inner field. The solid lines 
represent the LMC P-L relations derived by Udalski et al. (1999), shifted to the appropriate mean relative distance modulus. The dashed 
lines indicate the 2<t dispersion. 
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Fig. 15. — Relative distance modulus versus period for the restricted sample of Cepheids in the outer field. Crosses: Cepheid candidates 
that were rejected by selection criteria 1-4. Open circles: Candidates rejected by period cut or distance modulus clipping. Filled circles: 
Final sample of Cepheids used to determine the mean relative distance modulus. Solid line: mean relative distance modulus. Dotted lines: 
la dispersion of the final sample. Dashed line: Final adopted period cut. A typical error bar is shown on one of the data points. 
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Fig. 16. — Relative distance modulus versus period for the restricted sample of Cepheids in the inner field. Crosses: Cepheid candidates 
that were rejected by selection criteria 1-4. Open circles: Candidates rejected by period cut or distance modulus clipping. Filled circles: 
Final sample of Cepheids used to determine the mean relative distance modulus. Solid line: mean relative distance modulus. Dotted lines: 
la dispersion of the final sample. Dashed line: Final adopted period cut. A typical error bar is shown on one of the data points. 
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Fig. 17. — Mean distance modulus (relative to the LMC) as a function of cut at minimum period, for the Ly > 2 samples of the outer 
(top) and inner (bottom) fields. The error bars represent the la uncertainty in the mean. Our final choices for minimum period cut are 
indicated with open circles. 
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Fig. 19. — Residual of the individual Cepheid distance moduli about the best-fit line of Figure 18, plotted as a function of E{B — V). 
Cephcids located in the outer field are represented by concentric open and filled symbols. Cepheids located in the inner field are indicated 
by filled symbols. 
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Fig. 21. — Cepheid metallicity dependence determined through a comparison of TRGB and Cepheid distance moduli for 20 fields in 18 
galaxies. This Figure reproduces the bottom panel of Figure 12 of SaJcai et al. (2004), with the addition of two points for our fields in 
NGC 4258, shown with concentric open and filled symbols. 
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Fig. 22. — Correlation between period and extinction-corrected LMC-relative distance moduli for different choices of input P-L relation. 
Top panel: Outer field Ccplicids. Bottom three panels: Inner field Cepheids. U99=Udalski ct al. (1999); S04=Sandage et al. (2004); 
T03=Tammann et al. (2003). The LMC P-L relations are good fits to the Cepheid samples of both fields. The adoption of the Milky-Way 
P-L relations of Tammann et al. (2003) leads to a residual slope with a significance of ~ 7<t level for Pmin = 6 d or ~ 4cr for Pmin = 12 d. 
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Fig. 23. — WMAP 3-year 1 and 2<t error contours (dashed lines) in the Qm — w plane, for the wcdm+nopert model of Spergel et al. (2006). 
The solid contours represent the improvement obtained by using priors on Ho. Left panel: prior of Hq = 72±7 km Mpc~^ (Preedman 
et al. 2001). Right panel; prior of hypothetical future measurement of Hq = 74 ± 3.5 km Mpc~^. 
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Fig. 24. — Icr uncertainty in the value of w for the wcdm+nopert model of Spergel et al. (2006), after including hypothetical priors on Hq 
with a value of 74 km s~^ Mpc~^ and decreasing uncertainty. 
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Table 1 
Log of Observations 



Visit UT Date MJD at mid exposure* 

# B I V B I V 



O-Ol 2003 December 6 2980.0640 0.0734 0.0823 0.1206 0.1300 0.1389 
O-02 2003 December 7 2980.7980 0.8073 0.8163 0.8533 0.8627 0.8716 
Table appears in its entirety in the full-resolution version of the paper 

Note. *: JD-2450000.0 for first exposure; thereafter, only the last five digits are given, f : Guide star problems; 
limited usefulness. O: outer field; 1: inner field. 



Table 2 
Secondary standards 



ID R.A. Dec. X Y V I B 

(J2000.0) (pix) (mag) 

0-42521 12:19:12.170 47:10:55.12 3587.1 4181.2 23.203(10) 23.099(16) 23.275(05) 
0-42199 12:19:12.218 47:10:48.86 3704.5 4137.1 23.171(07) 22.911(07) 23.367(06) 

Tabic appeal's ui iJ.s cnlircly iii Uic full-rcHolulion dcthioii of Uic paper 

Note. — Uncertainties are given in parentheses and are expressed in units of 10~^ mag. O: outer field; I: inner 
field. 



Table 3 
Effects of selection criteria 



Selection criteria (§3.3) 



Initial Sample 

1. Amplitude ratios 

2. E{B-V) < foreground (2a) 

3. E{B-V) > 0.5 mag 

4. A/zo > 12 mag & a clipping 

Final sample 



Lv > 0.75 Lv > 2.0 



Inner 


Outer 


Inner 


Outer 


402 


134 


195 


63 


-110 


-28 


-23 


-4 


-37 


-30 


-17 


-20 


-9 


-1 


-6 


-1 


-28 


-12 


-11 


-3 


218 


63 


138 


35 



Note. — This table shows the reduction in size of the different samples as a result of the selection criteria listed 
in §3.3. 
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Table 4 

Cepheid variables - Basic data 



ID Per. R.A. Dec. X Y 

(d) (J2000.0) (pixels) V 


Magnitudes 
I 


B 


Amplitudes 
V I B 


Lv 


0-15165 3.30 12:19:20.668 47:10:31.19 3567.3 2383.1 26.646(042) 
0-21998 3.36 12:19:17.808 47:10:11.81 4100.5 2836.1 26.540(043) 
Table appears in its entirety in the full-resolution version of the paper 


26.092(040) 
25.959(043) 


27.148(054) 
26.993(061) 


503 168 691 
482 214 812 


2.01 
2.26 



Note. — Errors in mean magnitudes arc shown in parenthesis and expressed in units of 10 ^ mag. Light curve 
semi-ampHtudes are expressed in units of 10~^ mag. O: outer field; I: inner field. 



Table 5 
Rejected Cepheid candidates 



id Per. R.A. Dec. X Y 

(d) (,72000.0) (pixels) V 


Magnitudes 
I 


B 


Amplitudes 
V I B 


Lv 


Nv < 18 or Ni < 18 or JVb < 18 

0-18637 3.21 12:19:21.128 47:11:27.81 2453.4 2607.1 27.000(051) 
O-03121 3.23 12:19:27.280 47:10:20.74 3394.4 1030.0 27.158(030) 
Table appears in its entirety in the full-resolution version of the paper 


26.281(052) 
26.623(113) 


27.409(077) 
26.981(027) 


361 190 474 
359 29 


1.27 
0.92 



Note. — Errors in mean magnitudes are shown in parenthesis and expressed in units of 10 ^ mag. Light curve 
semi-amplitudes are expressed in units of 10~^ mag. O: outer field; I: inner field. 



Table 6 

Cepheid variables - Derived properties 



ID Per. E{V-I) /zg^ E{B-V) r/ 

(d) (mag) (mag) (mag) (mag) nso 

0-15165 3.30 11.009(115) 0.022(025) 11.009(075) 0.016(025) 1.30 
0-21998 3.36 10.900(122) 0.022(025) 10.901(075) 0.016(025) 1.23 
Table appears in its entirety in the full-resolution version of the paper 

Note. — l^o^'- Extinction corrected distance modulus derived from VSzI data (Eq. 7). /ig^: Extinction corrected 
distance modulus derived from the average of Eqns. (7-9). E{B — V): Average value of extinction derived from all 
measured color excesses. Errors in distance moduli and extinction are shown in parenthesis and expressed in units of 
10~^ mag. O: outer field; I: inner field. 
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Table 7 

Error Budget of the Cepheid Distance Scale 





1— ^Vi^TT"! rAl 1 C 

JT i L- V i\J Llo 


This work 


Goal 


A. Fiducial galaxy 


LMC 


NGC4258 


SI. Distance modulus (sys) 


0.13 


0.12 


0.04 


Rl. Distance modulus (ran) 




0.09 


0.02 


B. Photometric calibration 








S2a. V zcropoint 


0.03 


0.02 


0.02 


S2b. / zeropoint 


0.03 


0.02 


0.02 


S2. Photometry (sys) 


0.09 


0.05 


0.05 


R2. Photometry (ran) 


0.05 


0.03 


0.02 


C. Extinction corrections 








R3. Uncertainty in Ry 


0.02 


0.02 


0.02 


R4. Dc-rcddened PL fit 


0.04 


0.02 


0.02 


D. Metallicity corrections 








S3. Adopted correction 


0.08 


0.04 


0.03 


Rt- Total random 


0.07 


0.10 


0.04 


St- Total systematic 


0.18 


0.14 


0.07 


Combined error (mag) 


0.19 


0.17 


0.08 


Combined error (%) 


10 


8 


4 



Note. — All errors expressed in magnitudes unless otherwise indicated. Previous: adapted from Gibson et al. 
(2000). Goal: Anticipated reduction in uncertainties from Humphreys et al., (in prep.), Bersier et al., (in prep.) 
and Macri et al., (in prep.). 



Table 8 

Updated distance moduli to high-quality type Ia SNe 



Galaxy 


SN 


[0/H] 






MO 


Ref. 


name 




(dex) 


(mag) 


(mag) 


(mag) 




NGC 3370 


1994ae 


8.80 ±0.05 


32.23 ±0.04 


32.31 ±0.06 - 


-19.15 ±0.12 


R05 


NGC3982 


1998aq 


8.75 ± 0.05 


31.56 ±0.08 


31.63 ±0.09 - 


-19.15 ±0.12 


SOI 


NGC 4536 


1981B 


8.85 ±0.20 


30.80 ±0.04 


30.90 ±0.06 - 


-19.18 ±0.12 


FOl 


NGC 4639 


1990N 


9.00 ± 0.20 


31.61 ±0.08 


31.75 ±0.09 - 


-19.08 ±0.12 


FOl 


Average 










= -19.14 ±0.07 



Note. — ^o- Published Cepheid distance moduli; /xo,z : Distance moduli corrected for metallicity and our 
determination of the distance to the LMC. References: FOl = Freedman et al. (2001); R05 = Riess et al. (2005); SOI 
= Stetson & Gibson (2001). 
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Table 9 
Cepheid photometry 



Visit V I B V I B V I B 



0-15165 P = 3.30 d 0-21998 P = 3.36 d O-10450 P = 3.71 d 

1 26.900(096) 26.247(165) 27.532(121) 26.878(183) 25.915(120) 27.575(117) 26.565(116) 25.775(154) 26.898(109) 

2 26.797(163) 26.644(287) 26.448(325) 27.048(150) 26.091(208) ... 26.749(151) 26.068(123) 27.192(118) 

3 27.160(118) 26.413(162) 27.715(199) 27.090(160) 26.795(260) 27.453(272) 26.981(232) 26.310(219) 27.634(181) 

4 27.000(188) 26.367(161) 27.304(157) 27.178(207) 26.153(187) 27.661(365) 26.955(160) 26.114(163) 27.653(244) 

5 26.118(073) 25.837(117) 26.391(076) 25.992(062) 25.825(114) 26.137(068) 26.912(091) 26.472(214) 27.912(236) 

6 26.121(084) 25.973(144) 26.370(146) 25.942(089) 25.873(229) 26.199(104) 27.345(194) 25.990(117) 27.864(277) 

7 27.430(268) ... 27.672(170) 26.683(115) 26.082(126) 27.116(139) 26.127(088) 25.901(191) 26.461(093) 

8 26.994(198) 25.925(142) 27.453(167) 26.887(169) 26.252(180) 27.225(180) 26.268(083) 25.710(155) 26.591(144) 

9 26.283(084) 25.768(083) 26.304(082) 26.033(098) 25.567(120) 26.242(072) 26.607(127) 26.262(250) 27.656(200) 

10 25.882(135) 25.868(165) 26.605(093) 26.074(097) 25.677(095) 26.222(064) 27.255(166) 25.975(155) 27.197(212) 

11 27.057(247) 26.516(307) 27.971(295) 27.290(340) 26.589(342) 28.269(270) 26.768(186) 25.892(163) 26.691(078) 

12 27.188(139) 26.129(103) 27.722(315) 26.773(105) 26.177(132) 27.532(207) 26.223(044) 26.346(206) 26.622(091) 

13 26.727(171) 26.342(221) 27.554(137) 26.749(111) 25.874(165) 27.785(287) 27.022(256) 26.507(205) 

14 26.691(188) 25.951(117) 28.008(229) 26.812(142) 26.442(275) 27.492(234) 26.349(072) 25.814(135) 26.683(122) 



Note. — The Julian Date for each visit can be found in Table 1. (*): Measurement deviated by more than 
from best-fit template light curve and was rejected. This table is available in its entirety upon request. 



